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Abstract. The ultra-relativistic heavy-ion program at RHIC has shown that at intermediate trans- 
verse momenta (pr — 2-6 GeV) standard (independent) parton fragmentation can neither describe 
the observed baryon-to-meson ratios nor the empirical scaling of the hadronic elliptic flow (V2) ac- 
cording to the number of valence quarks. Both aspects find instead a natural explanation in a coales- 
cence plus fragmentation approach to hadronization. After a brief review of the main results for light 
quarks, we focus on heavy quarks showing that a combined fragmentation and quark-coalescence 
framework is relevant also here. Moreover, within relativistic Langevin simulations we find evi- 
dence for the importance of heavy-light resonances in the Quark-Gluon Plasma (QGP) to explain 
the strong energy loss and collective flow of heavy-quark spectra as inferred from non-photonic 
electron observables. Such heavy-light resonances can pave the way to a unified understanding of 
the microscopic structure of the QGP and its subsequent hadronization by coalescence. 
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INTRODUCTION 

Experimental data from the Relativistic Heavy Ion Collider (RHIC) during the past few 
years have shown convincing evidence for a state of matter with energy densities, £, in 
substantial excess of the expected critical one, i.e., £ ~ 15-20e c with £ c ~ 1 GeV/fm -3 . 
In the standard picture before 2003 the main concern in the study of produced hadrons 
at high transverse momentum (pr) was the underlying modification of partonic spec- 
tra due to interactions in such a hot and dense medium, that would reflect in a similar 
pattern for all hadrons through independent fragmentation. Indeed, one of the most ex- 
citing observations at RHIC has been the suppression of high-p^ particles in agreement 
with the non-abelian radiative energy-loss theory within perturbative QCD (pQCD) |H]. 
However, for light hadrons the observations of an anomalous baryon-to-meson produc- 
tion ratio at intermediate pj up to ~ 6 GeV and a scaling of the elliptic flow with the 
number of quark constituents, has enforced revisions of an independent fragmentation 
model for hadronization. Instead coalescence processes among massive quarks appear 
to be a convenient picture that can naturally and quantitatively account for the main 
features of light hadron production at intermediate pr iSBH,^. 

For heavy quarks (charm (c) and bottom (b)), the energy loss predicted by pQCD [@] 
turned out to be insufficient (at variance with the light quark case) to account for 
the observed large nuclear suppression (small Raa) and collectivity (large V2) in non- 
photonic single-electron spectra [Q, 0, S Gil [Hfl. Here the challenge is mainly in the 



understanding of the in-medium quark interactions, even if the acquired knowledge on 
the hadronization mechanism from light quarks plays a significant role as well. 

Lattice QCD (1QCD) results suggest that resonance structures survive in meson- 
correlation functions at moderate temperatures 112D above T c . Two of us have there- 
fore suggested an effective model for heavy-light quark scattering via D and B res- 
onances |131]. To study the consequences of such interactions for the modification 
of the heavy-quark (HQ) distributions, a Langevin simulation has been performed to 
trace the evolution of HQ distributions through the fireball in heavy-ion Collisions 
(HICs) [EUdl; hadronization has been modeled by a coalescence model similar to the 
one applied to light quarks [0,1 . It has been found that the effect of resonant heavy-light 
scattering is crucial QJL5h and also provides a reasonable agreement with semileptonic e ± 
spectra at RHIC dataB O, Gil • 

We note that for heavy quarks (heavy-light) quark-antiquark resonances provide the 
dominant medium effects on their distributions close to T c , which then naturally merges 
into a coalescence-type description for hadronization processes jlol . 



MODIFICATION OF HADRONIZATION MECHANISM 

Hadronization at asymptotically large momentum can be described by a set of fragmen- 
tation functions D a / H (p/P) that parametrize, in a universal way, the probability that a 
hadron H with momentum P is created from a parton a with momentum p in the vac- 
uum. Fragmentation functions have been measured in e + e collisions and work well for 
hadron production at pj> 2 GeV also in pp collisions at RHIC energies. Therefore it was 
expected that in this pj regime the QGP could be probed by focusing on modifications 
of the spectra, EdN/d 3 p, at the parton level; but from Au+Au collisions it became clear 
that this is not the case at least up to pj ~ 6 GeV. Two puzzling observations lead to this 
conclusion: (a) baryons are much more abundant than predicted by fragmentation. A ra- 
tio pj% ps 1 between 2 and 4 GeV/c has been measured, see Fig.[H much larger than the 
value of ps 0.2 predicted by leading-twist pQCD. A similar trend is observed for p/n, 
A/K^, i.e., p, p and A's do not seem to suffer jet quenching. The pertinent nuclear mod- 
ification factors, Raa, are close to 1, unlike those of light mesons for which Raa ~ 0.2 
in central collisions; (b) the elliptic flow of all identified hadrons is found to scale ac- 
cording to a quark-number scaling as reflected in a universal behavior of V2h(pt / n )/ n 
where n is the number of constituent quarks in hadron H. In particular, recent data for the 
0(1020) also follow the scaling, suggesting the dominance of the quark content rather 
than the mass effect, in agreement with the coalescence prescription [17]. 

The main reason for the inadequacy of a pure fragmentation picture is the high 
density of the matter created in HICs. In such an environment one may expect that 
quarks could just coalesce into hadrons: three quarks into a baryon, a quark-antiquark 
pair into a meson. In such a picture baryons with momentum pr are mainly produced 
from quarks with momenta ~p T /3, while mesons with the same momentum mainly arise 
from quarks with momenta ~p T /2. This is contrary to the fragmentation process where 
baryon production is suppressed with respect to mesons as more quarks are needed from 
the vacuum. A coalescence model that is based on the simple overlap of the quark- 
distribution function with a hadron- wave function has been developed to implement the 
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FIGURE 1. Left: it and p~ transverse-momentum spectra in v / i=200 AGeV Au+Au collisions. RHIC 
data Ill8lll9ll are shown by circles (n°) and squares (p), lines are results from coalescence. Right: experi- 
ment vs. coalescence results for baryon-to-meson ratios for p/n (lower part) and A/K® (upper part). 



physical ideas sketched above [Q]. In such a model the transverse-momentum spectrum 
of hadrons that consists of n (anti-) quarks is given by the overlap between the hadron- 
wave function and n quark phase-space distribution functions, f q (xi,pi): 



~£~p^ =gH j EI (27l)lE P ' - A °ifq^Pi)ffi{ x \-- x n\Pi--Pn) 5 (2) \ ~Y,PT,i) \ 



(1) 



dcr denotes an element of a space-like hadronization hypersurface, fn is the Wigner 
distribution function of the hadron and gjj the probability of forming a color-neutral 
object with the spin of the hadron from n colored quarks. Therefore it is assumed that 
the probability of coalescence is simply given by the phase-space distance weighted by 
the wave function of the produced particle Eq. (Q~|). In such an approach constituent- 
quark masses are included representative for non-perturbative effects. This is a further 
assumption that can be relaxed for heavy quarks owing to the smaller effect of the 
QCD vacuum on their masses. The distributions f q (xi,qi) are fixed as homogeneous 
Boltzmann distributions with an average radial flow (j8) =0.35 for pj < 2 GeV, and 
quenched minijets for pj > 2 GeV. The volume is fixed to reproduce the measured 
transverse energy at given centrality. In the left panel of Fig. [Q the p and % spectra 
obtained from the quark-coalescence model are shown together with the experimental 
data from PHENIX [H . 



The resulting p/n ratio is shown in the right panel of Fig. Q] by lower lines D20Q, a 
similar effect is seen also for the A/K® ratio (upper lines) 1121 1 . Similar conclusions are 
reached in other studies based on quark coalescence 0,0]. 



HEAVY-QUARKS AT HIGH TEMPERATURES 

Heavy quarks (£>,c) are produced out of thermal equilibrium in the very early stage 
of the reaction; due to their large mass a perturbative evaluation of their in-medium 
interactions was expected to be reliable also at relatively small pj. However, a small 
nuclear modification factor, Raa — 0.3, has been deduced from semileptonic electron 



spectra associated with decays of D- and 5-mesons [81, HID, comparable to the pion 
one. Such a value is incompatible with pQCD jet-quenching mechanisms [6i]. This 
statement is strengthened by the observed V2 of up to 10% nLLLOO, indicating substantial 
collective behavior of charm (c) quarks yfl. Moreover, a consistent description of Raa 
and V2 cannot be achieved even if one artificially upscales the transport coefficients 
within pQCD energy-loss calculations jl4 \. This suggests that the physics underlying 
the heavy quark observables is not only a matter of a global evaluation of the interaction 
strength, but there is an opportunity for a more detailed understanding of the microscopic 
nature of the interaction (most likely of non-perturbative origin) and of the hadronization 
mechanism, as we briefly review in the following. 

A hint on non-perturbative interactions of heavy quarks in the medium is provided by 
1QCD computations which exhibit resonance structures in meson correlation functions 



at moderate temperatures | |12L l22 |l . Along this line two of us have suggested that D- and 



5-resonance exchange in the q — Q channel may be the dominant scattering process that 



drives the HQ dynamics [13]. To evaluate the consequences of such a picture we have 



built a model based on an effective Lagrangian: 

^ = Q^-^Tq + \v.c. (2) 

where <&=D, B. We have calculated elastic Q + q — > Q + q scattering amplitudes via 4> 
exchange in the s- and w-channel. The existence of one ^> state (e.g., a pseudoscalar 
J p =0~), is assumed together with a minimal degeneracy following from chiral and HQ 
symmetries, represented by Dirac matrices T=l, 75, 7^, 757^ in Eq. ©. 

The application to HICs is realized by treating HQ kinetics in the QGP as a relativistic 



Langevin process [1511 



df _ djypf) | d 2 (D p f) ; (3) 
dt dp dp 2 

7 and Dp are drag and (momentum) diffusion coefficients which determine the ap- 
proach to equilibrium and satisfy the Einstein relation, T = D p / jMq. The bulk medium 
is modeled by a spatially homogeneous elliptic thermal fireball which expands isen- 
tropically. Finally, hadronization is treated via coalescence at r c =180 MeV, see Eq. (0Q), 
plus fragmentation processes evaluated as / C) &(pr) * [1 — -P c ,fc->(D,A c ),(£,A&)(/ 7 r)]> where 
P c ,b^(D,A c ),(B,A b ) is m e probability for a heavy quark to coalesce. 

For 200 AGeV Au+Au collisions results from the Langevin simulation including 
hadronization by coalescence+fragmentation (left) and fragmentation only (right) are 
shown in Fig. [Altogether with experimental data [8j,|9|]. It is obvious that elastic scattering 
in a pQCD scheme is insufficient to account for the small Raa, independent of the 
hadronization scheme applied. The red band shows the full calculation with c, b quarks 
that scatter in the presence of hadron-like resonances with a width T ~ 0.4-0.75 GeV 
(representing the interaction strength), supplemented by the pQCD elastic scattering in 
color non- singlet channels (dominated by gluons). We note that the contamination of 
single electrons from B decays is significant already at pj ~ 2 GeV (corresponding 
to a crossing of c and b spectra at around pj~4-5 GeV). Thus the inclusion of B 
mesons (despite the inherent uncertainties in the blc ratio) is mandatory to draw reliable 
conclusions on the interaction processes underlying the experimental results. 
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FIGURE 2. Nuclear modification factor for single electrons, including coalescence and fragmentation 
at hadronization (left panel) and only with fragmentation (right panel), compared to RHIC data lflo[[Tlll . 
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FIGURE 3. Elliptic flow for single electrons, including coalescence and fragmentation at hadronization 
(left panel) and only with fragmentation (right panel), compared to RHIC data Uoll . 



When comparing the band in the left and right panel of Fig. [2] we can see a clear 
coalescence effect of hadronization in terms of an increase in the Raa at Pt ~ 1-4 GeV. 
This effect becomes more significant via a simultaneous enhancement of the elliptic flow 
V2 [0, [3], see Fig. [3] This behavior is typical of a coalescence mechanism, reversing 
the usual correlation between Raa and vi, and allows for a reasonable description of the 
experimental data on both Raa and V2- Very recently, a Brueckner many-body scheme for 
in-medium T-matrices for HQ scattering off light quarks [12311 has been evaluated starting 
from 1QCD potentials. The existence of D and 5-like resonances has been confirmed; 
furthermore, when embedding the pertinent drag and diffusion coefficients in the model 
described above, a comparable (or even better) agreement with the data is found [23]. 



CONCLUSIONS 

The first stage of RHIC program has shown clear signs of modifications of the hadroniza- 
tion mechanism in the light-quark sector relative to pp collisions. There are several ev- 
idences that hadronization proceeds through the coalescence of (massive) anti-/quarks 
which are close in phase space. The modification of hadronization for light quarks also 



seems to play a role in the new challenge posed by heavy-quark probes. Here, the main 
issue is the dominant interaction mechanism (HQ diffusion) and its relation to the mi- 
croscopic structure of the QGP. We have presented a scenario based on the existence of 
B- and D-like resonances in the QGP up to T = 2 T c . Reasonable agreement with exper- 
imental data is achieved owing to a positive interplay between the resonance scattering 
mechanism and the hadronization by coalescence. We note that such an approach pro- 
vides an inherent consistency between the in-medium interactions of HQs and the sub- 
sequent hadronization by coalescence: a pole in the quark-antiquark propagator above 
T c can be viewed as a precursor of recombination. 

Finally, it is important to keep in mind the interrelations of open and hidden charm, 
if regeneration makes up a good fraction of the charmonium yield as expected at RHIC 
and LHC. In this case, the study of B and D distributions can be related to that of 7/ i// 
and T due to the underlying common c, b distributions; the consistency of such a picture 
should be checked in the near future. 
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